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ABSTRACT 

Characterization of exoplanets has matured in recent years, particularly through studies of exo- 
planetary atmospheres of transiting planets at infra-red wavelenegths. The primary source for such 
observations has been the Spitzer Space Telescope but these studies are anticipated to continue with 
the James Webb Space Telescope (JWST). A relatively unexplored region of exoplanet parameter 
space is the thermal detection of long-period eccentric planets during periastron passage. Here we 
describe the thermal properties and albedos of long-period giant planets along with the eccentrici- 
ties of those orbits which allow them to remain within the habitable zone. We further apply these 
results to the known exoplanets by calculating temperatures and flux ratios for the IRAC passbands 
occupied by warm Spitzer, considering both low and high thermal redistribution efficiencies from the 
perspective of an observer. We conclude with recommendations on which targets are best suited for 
follow-up observations. 

Subject headings: planetary systems - techniques: photometric 



1. INTRODUCTION 

As the number of known exoplanets continues to rise 
at a steady pace, their diversity appears to only increase. 
For instance, transiting exoplanets have allowed us to ac- 
cess properties of exoplanet atmospheres through obser- 
vations at Infra-Red (IR) wavelengths during secondary 
ecli pse and atmospheric absorption during primary tran- 
sit (jAgol et alj|2010t iDeming et al.ll2007at) . Further op- 
portunities for atmospheric studies have arisen through 
the detection of phase variations of such planets as 
HP 189733b (iKnutson et al.1 l2009ai) and HD 149026b 
(|Knutson et all l2009bl ). These phase variations have 
also been detected for non-transiting planets, includin, 



v And b (ICrossfield et al.ll2010t lHarrington et alll20(> 
and HD 179949b (|Cowan et all 120071). Significant con 



straints have been placed on the planet-to-star flux ratio 
for HD 217107b through grou nd-based near-IR observa- 
tions bv lCubillos eTaflll20ll . 

Studies thus far have largely been directed towards 
short-period planets which are expected to have an a 
priori high effective temperature. In particular, this ex- 
clusivity results from a relatively high thermal flux re- 
quired from the planet in order for the signature to be de- 
tectable. The exception to this is the planet HD 80606b, 
with a highly eccentric o rbit (e = . 93) an d a period of 
111 days, detected bv lNaef et all (120011). Su bsequent 



Spitzer observations by iLaughlin et all ()2009l ). as well 
as the detection of the secondary eclipse of the planet, 
were used to measure the out-of-eclipse variations and 
estimate the ra diative time constant at 8 microns. As 
pointed out by iBarbieri et al.1 ()2007f ) , ambiguous mea- 
surements of radiative time constants have prevented 
a consensus on expected planet-wide flow patterns for 
short-period planets. 

Recent three-dimensional models of planetary atmo- 
spheres dBarman et ail 120051: iDobbs-Dixon et all 120101: 
iFortnev et al.1 120101; iKoskinen et all 120071 iLewis et al. 
120101: [Madhusudhan & Scagcr 2010; Rausch er fc Menou 



20101: IShowman et all [20081: ISudarskv et all 120051: 
Thrastarson fc Choi 120101) " have made significant 
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progress in deriving the underlying atmospheric physics 
which drive the thermal properties and zonal winds in 
exoplanetary atmospheres. Highly eccentric planets such 
as HD 80606b provide a means to explore atmospheric 
properties in a different regime of orbital parameter 
space since the heating of the atmosphere during peri- 
astron passage can be sufficien t to produce a detectable 
sis nal ([C owan fc Agoll I2011al ). The possability of 
investigating eccentric planets at optical wave l engths 
has previously been explo red by iKane fc Gelinol (|2010f ) 
and IKane fc Gelinol (|2011l) . This observing window for 
periastron passage is brief and requires a reasonable 
unde rstanding of the orbital parameters (jKane et all 
I2009T) . 

This paper describes the predicted thermal changes for 
exoplanets in highly eccentric orbits. This study is pri- 
marily motivated from an observers point of view and 
is mostly concerned with planets which are not known 
to produce either a secondary eclipse or primary tran- 
sit. We derive analytical expressions for the albedos 
based upon theoretical models and calculate the effec- 
tive temperatures and flux ratios, taking into account the 
thermal heat redistribution and radiative time constant. 
These calculations are applied to the known exoplanets 
for which orbital parameters measured from radial veloc- 
ity data are available. We explore the dependencies of 
the planetary effective temperatures on eccentricity and 
orbital period and determine the percentage of the or- 
bits which are spent in their respective habitable zones. 
We finally calculate predicted maximum flux ratios dur- 
ing periastron passage, discuss the effect of spots on de- 
tections, and propose potential targets for Spitzer and 
JWST observations. 

2. PLANETARY EFFECTIVE TEMPERATURES 

In this section we outline the basic assumptions used to 
calculate the planetary effective temperatures through- 
out the remainder of the paper. These assumptions have 
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been deliberately introduced to be fairly broad since the 
intention is to encompass a variety of planets to pro- 
duce a first-order approximation of the global distribu- 
tion of thermal signatures and detectability. These were 
designed in such a way so as to produce testable limits 
on the flux ratios at periastron passage. 

We begin with the luminosity of the host star, which 
is approximated as 

L* = 47ri^aT e 4 ff (1) 

where R± is the stellar radius, T e s is the stellar effective 
temperature, and a is the Stefan-Boltzmann constant. In 
cases where the radius of the star is not available from 
direct measurements, we estimate the radius from the 
derived values of the surface gravity log g using the rela- 
tion 



+ log 5© (2) 



where logp = 4.4 374 (ISmallevil2005l). 

As described by iKnutson et al.l (|200 9b) , we can ap- 
proximate the effective temperature of a planet, T p , as 
a blackbody. This approximation will deviate slightly 
from the true temperature depending upon albedo, at- 
mospheric properties, and internal heating. Assuming 
that the atmosphere is f00% efficient at redistributing 
heat around the planet, the planetary equilibrium effec- 
tive temperature is given by 



X = 

lQirar 2 



(3) 



where A is the spherical (Bond) albedo and r is the 
star-planet separation. In this case the surface is uni- 
formly bright and thus there will be no observable phase 
function at infra-red wavelengths. However, if the atmo- 
sphere is inefficient with respect to heat redistribution, 
this will lead to a hot dayside for the planet where the 
effective temperature is 



L*(l - A) 
8nar 2 



(4) 



where there will be a resulting phase variation as the 
planet orbits the star. The generalized form for the plan- 
etary effective temperature is thus 

_ / L^l-A) V 
p ~ ^(l + 77)87rar 2 ) [ ) 

where r\ is the atmospheric heat redistribution efficiency 
with a value ranging between and f . 

For a typical hot- Jupiter scenario, the star-planet sep- 
aration is assumed to be the same as the semi-major 
axis, a, since these are usually circular orbits. However, 
the star-planet separation for eccentric planets has the 
following form: 

r= a(1 - e2) , (6) 
1 + ecos/ w 

where e is the orbital eccentricity and / is the true 
anomaly. Thus, the eccentricity of a planetary orbit in- 
troduces a time dependency to the effective temperature 
of the planet. In the absence of atmospheric effects (see 
Section [3]) , the temperature of the planet may be exam- 
ined as indicated in Figure [I] where the dependence on 
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Fig. 1. — Dependence of planetary effective temperature on star- 
planet separation for inefficient (rj = 0) and efficient (rj = f) heat 
redistributions. This is calculated for GOV and K5V host stars and 
shows that the scale of the temperature variations are dominated 
by the star— planet separation and not by the heat redistribution 
efficiency. 

star-planet separation is shown for two example spectral 
class targets. For a given spectral type, the difference 
between 0% and 100% heat redistribution efficiency can 
lead to a ~ 20% adjustment in the temperature calcula- 
tion. However, r is clearly dominant over r\ at all sep- 
arations. This means that the planetary temperature, 
and thus the flux ratio, is more dependent upon the or- 
bital properties (which are well-determined) than on the 
atmospheric dynamics (which are generally unknown). 
This is elaborated upon in the following section. 

3. ATMOSPHERIC PROPERTIES 

Here we discuss the major atmospheric properties 
which directly influence planetary thermal signatures. 

3.1. Spherical Albedo 

The spherical albedo of planetary atmospheres be- 
yond the hot Jupiter regime is less understood at IR 
wavelengths (~ 10 /im) than it is at optical w ave- 
lengths. However, the models of lMarlev et al.1 (|1999f ) and 
ISudarskv et al.l (|2005l ) indicate that the albedo drops sig- 
nificantly in the range between 1.0 /iin and 1.2 /im in 
a manner which is relatively independent of the star- 
planet separation. In the case of Jupiter, the effective 
temperature of Jupiter beyond ~ 5 /jm is ~ 125 K and 
the spherical albedo is ~ 60% of the geometric albedo. 
In addition, constraints p laced upon the s pheric al albe- 
dos of hot Jupiters (e.g., IKnutson et al.l ()2009aD ) show 
these to be exceptionally low and in agreement with mod- 
els which predict the removal of reflective cloud layers 
in those extreme stellar flux regimes. We thus gener- 
alize the al bedo dependenc e on s tar-planet separation 
adopted by iKane fc Gelinol (|2010l) by scaling the rela- 
tion for spherical albedos. However, here we are testing 
the conditions at periastron where the albedo will be 
minimum. This is justified by the low measured albedos 
mentione d above and th e indep endent statistical verifi- 
cation bv lCowan fc Agol ()2011b| ) which favors low Bond 
albedos for small star-planet separations. 

3.2. Heat Redistribution Efficiency 

The small phas e ampli tude of HD 189733b detected 
by IKnutson et al.l (|2009al) indicates that this particular 
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planet has very high heat redistribution efficiency caused 
by atmospheric advection which produces high-speed 
zonal winds which carry heat to the night-side of the 
planet. These atmospheric pa tterns are highly model de- 
pendent ([Barbieri et al.ll2"007l ) and poorly understood for 
planets in the short-period regime due to the small sam- 
ple size and the complex interaction of planetary struc- 
ture, composition, tidal effects, and incident flux. Longer 
period eccentric orbits such as HD 80606b will further 
exhibit time-dependent behaviour depending upon the 
star-planet separation and ps eudo- s ynchr onized spin ro- 
tation (see equation (42) of iHud (j 19811 )). A statisti- 
cal study of 24 known transiting planets performed by 
ICowan fe Agoll (|2011b[ ) found that there is expected to 
be a wide range in heat redistribution efficiencies. As 
described earlier, the star-planet separation is domi- 
nant over heat redistribution efficiency in determining 
the planetary effective temperature. However, given the 
uncertainties in planetary models, we consider the two 
extremes of r\ — and n = 1 to determine the upper and 
lower bounds on the flux ratio for a given planet. 

3.3. Radiative and Advective Time Constants 

The radiative time constant is a quantitative mea- 
sure of the seasonal lag caused by t he thermal response 
of the atmosphere to incident flux (jFortnev et al.l [20081 : 
ISeager et al.ll2005l) . The radiative time constant, r ra d 
is related to fundamental atmospheric properties in the 
following way 

P c P 

Trad ~.g4aT3 { ' } 

where P is the pressure, g is the surface gravity, and 
cp is the specific heat capacity. If r ra( j = then the 
incident flux is immediately re-radiat ed from the day-side 
of the planet (|Cowan fe Ago] 1201 laf). For example, the 
measured value for HD 80606b bv lLaughlin et atl (|2009f ) 
is T ra d = 4.5 ± 2 hours. 

The related quantity is the advective time constant, 
Tadv , which is a measure of the movement of a parcel of 
gas around the planet. This is given by 



7"adv 



Rp 
U 



(8) 



where U is the characteristic wind speed. Thus, r a( jv ap- 
proximates to zero when the wind speed becomes large, a 
situation which results in a high heat redistribution effi- 
ciency, as discussed in Section [3~2l The predicted values 
of r ra d and r a dv are highly depen dent upon the circu- 
lation models of the atmos p heres (lLangton fe La ughlin 
120081 iMontalto et~aT] l201lL iShowman et al.l 120091 ) and 
also vary with the wavelength since this e ffects the depth 
into w hich the atmosphere is probed (Knutson et al. 
2009a). Here we concern ourselves with the peak flux 
ratio which is expected to occur relatively close to the 
point of periastron passage for eccentric orbits. Thus 
we assume that r ra d <C P and r a d v 3> P for the subse- 
quent calculations, keeping in mind that significant plan- 
etary spin, in addition to atmospheric composition and 
cloud effects, may induce wind patterns that cause di- 
vergence from this assumption. This produces an upper 
limit on the flux from the planetary substellar point and 
thus an upper limit on the predicted flux ratio with the 
star. The interplay between the radiative/advective time 



constants and other effects, such as drag mechanisms 
and numerical dissi pation , have been investig ated by 
iRauscher fe Menoul pOllh : iThrastarson fe Cnol(|2011ft . 

4. FLUX RATIO AT PERIASTRON PASSAGE 

The measurable quantity from observations acquired 
at frequency v is the flux ratio between the star and the 
planet, given by 



El 



(exp {hv/kT cS ) - l)R 2 p 
(exp {hv/kT p ) - l)Rl 



(9) 



where R p is the radius of the planet. The flux ratio as 
a function of the orbital phase depends upon the ther- 
mal redistribution efficiency of the atmosphere. If this is 
100% efficient then the flux depends purely on the star- 
planet separation since th ere is no longer a phase function 
such as that described by iKane fe Gelinol ((20101) . 

Observations of HD 189733b and HD 80606b indi- 
cate that t hese planets have relatively high redistr i butio n 
efficiency (jKnutson et al.l l2009at iLaughlin et~aU 12009ft . 
These planets are in very different kinds of orbits and 
so this high efficiency may be quite common amongst 
planets which experience either constant or intermittent 
periods of high stellar flux. The lack of atmospheric 
measurements, and thus models to explain observations, 
for eccentric planets further motivates the need for this 
study. An alternative hypothetical explanation is that 
hot Jupiters have high zonal winds and very little cloud 
layer, whereas eccentric planets are not irradiated in the 
same way and so retain some of their clouds, even during 
periast ron passage, thus redu cing the redistribution effi- 
ciency (jCowan fe Agolll2~011aj) . This could lead to phase 
variations in the thermal signature such that planets 
whose periastron argument is uj ~ 270° become the op- 
timal targets. Discriminating between these competing 
ideas requires further observations to resolve. The calcu- 
lations in the following sections thus represent testable 
assumptions which can constrain these postulates. 

5. HABITABILITY OF ECCENTRIC PLANETS 

The Habitable Zone (HZ) is defined as the range of cir- 
cumstellar distances from a star within which a planet 
could have liquid water on its surface, given a dense 
enough atmosphere. The various cr iteria for defin i ng the 
HZ has been described in detail bv lKasting et al.l (1993) 
an d further generalized as a fun ction of spectral typ e 
by lUnderwood et all ([2003 ;■ and Uones fe Sleepl (|2010D . 
In estimatin g the boundaries of the HZ, we utlize the 
equations of lUnderwood et al.1 ()2003f ) which relate the 
radii of the HZ inner and outer edges to the luminosity 
and effective temperature of the host star. Using the 
boundary conditions of runaway greenhouse and maxi- 
mum greenhouse eff ects at the inner and ou ter edges of 
the HZ respectively (jUnderwood et aLll2003|) . the stellar 
flux at these boundaries are given by 

S- mneI = 4.190 x 10" 8 T o 2 ff - 2.139 x 10" 4 T cff + 1.268 

S onteI = 6.190 x 10~ 9 T 2 ff - 1.319 x 10~ 5 T cff + 0.2341 

The inner and outer edgers of the HZ are then derived 
from the following 
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Habitable Zone 



Fig. 2. — Eccentric orbits (dotted lines) relative to the inner 
and outer boundaries of the habitable zone (solid lines). The low 
eccentricity planet is able to maintain a presence inside the habit- 
able zone, whilst the high eccentricity planet moves beyond both 
boundaries. 

where the radii are in units of AU and the stellar lumi- 
nosities are in solar units. 

The effect of orbital eccentrici ty on the habitabili ty of 
planets has been inve stigated by lAtobe et al.1 (|2004[ ) and 
iDressing et al.l (|2010f ) . Here we expand upon this by us- 
ing our effective temperature calculations to determine 
the amount of time a planet spends in the HZ. There are 
two types of habitable zone situations we consider: plan- 
ets in eccentric orbits that remain within the bounds of 
the habitable zone, and planets in eccentric orbits that 
cross the boundaries of the habitable zone. In the lat- 
ter case, the atmospheric relaxation timescales (governed 
by r ra d and r a(1 v, described in Section |33[) may moder- 
ate temperatures even during dramatic changes in star- 
planet separation. Note that long-period planets in circu- 
lar orbits are more likely to have a non-synchronous spin 
rotation and so T ra d and T a dv will have a much smaller 
correlation with orbital period. However, planets in ec- 
centric orbits will likely spend the majority of their or- 
bits within the HZ if apastron lies close to the outer edge, 
thus retaining the possability of pseudo-synchronous spin 
rotation. Note that it is poorly understood to what de- 
gree exoplanets in eccentric orbits will retain a pseudo- 
synchronous spin rotation within the HZ due to the many 
contributing factors. However, the evidence for these rel- 
ative influences ha ve been observed for m oons within our 
own Solar System (|Gladman et al.|[T996l) . Figure [2] shows 
two example eccentric orbits overlaid on a hypothetical 
HZ which depicts the two cases described above. In the 
following section, we apply these methods to determine 
the percentage of the orbit spent in the HZ for the known 
exoplanets. 

6. APPLICATION TO KNOWN EXOPLANETS 

Here we calculate predicted effective temperatures and 
flux ratios at periastron passage for a sample of the 
known exoplanets for which there are known orbital solu- 
tions. The orbital parameters of 390 planets, along with 



the host star properties, were extracted us ing the Ex- 
oplanet Data Explorer 1 (jWright et all 1201 1[ ). The data 
are current as of 17th January 2011. Planets for which 
the host star values of T c g and logg were not avail- 
able were excluded from the sample. The calculated 
flux ratios are highly dependent on the assumed radii 
of the exoplanets (see Equation^ The anomlous radii of 
short-period tr ansiting planets has been investigated by 
iLaughlin et~aTl ([20111 ) . However, iFortnev et all ([20071 ) 
have shown that, for a given planetary mass and compo- 
sition, planetary radii should not vary substantially be- 
tween orbital radii of 0.1-2.0 AU. However, different com- 
postions, particularly with regards to core versus coreless 
models, can lead to radii variations of ~ 25%. Since we 
are considering eccentric orbits of massive planets, we 
take the conservative approach of fixing the radius for 
each of the planets in this sample at one Jupiter radius, 
keeping the caveats mentioned above in mind. 

There are a variety of current and future space-based 
observatories that are capable of effectively monitoring 
known exoplanets for thermal signatures. The IRAC in- 
strument of the Spitzer Space Telescope had passbands 
centered at 3.6, 4.5, 5.8, and 8.0 microns 2 . For the James 
Webb Space Telescope (JWST), planned instruments in- 
clude MIRI which covers 5 to 27 microns 3 , and NIRCam 
which covers 0.6 to 5.0 microns 4 . We calculate predicted 
flux ratios using the 3.6 and 4.5 micron passbands of 
IRAC since the longer wavelengths are no longer avail- 
able during the warm phase of the mission. 

Table [T] shows the results of these calculations for the 
50 most eccentric planets in the sample. Included in the 
table are the orbital period, P, eccentricity, e, periastron 
argument, cj, the percentage of a full orbit spent in the 
HZ, tfjz, and the effective temperatures and IR flux ratios 
assuming both low (rj = 0) and high (77 = 1) heat redistri- 
bution efficiencies. We also assume an albedo of A = 0.0, 
meaning that the planet absorbs 100% of the incident 
flux (seelMadhusu dhan fc Seagerl (|2009f ) for a more thor- 
ough statistical analysis of these values) . This produces 
a higher effective temperature for the planet but is a 
reasonable assumption since models and measurements 
have shown that planets lose their reflective cloud layers 
at small star-planet separations (jKane fc Gelinol [20101 : 
ISudarskv et all 120051 ). For comparison, an albedo of 
A = 0.3 leads to a ~ 10% reduction in the planetary 
effective temperatures. Note that flux ratios for those 
planets at particularly long periods are approximately 
zero at both passbands. However, recall from Figure [T] 
that the flux ratios are mostly dependent upon the star- 
planet separation at periastron rather than rj. Thus, as 
one increases the orbital period and decreases the ec- 
centricity, the change in 77 becomes less important (also 
demonstrated later in the top- right panel of Figure [3]). 
The flux ratios calculated for HD 80606b are comparable 
to the 0.0010±0.0002 values measured bv lLaughlin et all 
(2009) at 8 /im and indicate that there is moderate heat 
redistribution efficiency of the atmosphere in this case. 
Some of the more interesting examples are discussed in 
detail in Section [51 

1 http://exoplanets.org/ 

2 http:/ /ssc. spitzer. caltcch.edu/irac/ 

3 http: / / www.stsci.edu/jwst/instruments /miri / 

4 http:/ /www. stsci.edu/jwst/instruments/nircam/ 
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TABLE 1 

Effective temperatures and IR flux ratios for known exoplanets. 
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2151.0 


0. 


641 


330.2 


11.9 


330.3 


0.0000443 


0.0003698 


277.8 


0.0000045 





0000592 


HD 131664 b 


1951.0 


0. 


638 


149.7 


11.1 


322.1 


0.0000384 


0.0003408 


270.9 


0.0000037 





0000521 


HD 74156 b 


52.0 


0. 


630 


174.0 


0.0 


1235.3 


0.2161417 


0.3191066 


1038.7 


0.1150297 





1 895808 


HD 44219 b 


472.3 


0. 


610 


147.4 


73.9 


552.5 


0.0045355 


0.0143141 


464.6 


0.0011534 





0047790 


HD 154672 b 


163.9 


0. 


610 


265.0 


0.0 


766.4 


0.0358256 


0.0760486 


644.4 


0.0133102 





0342433 


HD 16175 b 


990.0 


0. 


600 


222.0 


23.8 


419.6 


0.0004756 


0.0023800 


352.8 


0.0000784 





0005626 


HD 3651 b 


62.2 


0. 


.596 


245.5 


0.0 


833.3 


0.1188941 


0.2310287 


700.7 


0.0477421 





1 1 05272 


HD 171028 b 


550.0 


0. 


590 


304.0 


54.7 


623.3 


0.0045460 


0.0121809 


524.1 


0.0013496 





0045978 


HIP 2247 b 


655.6 


0. 


540 


112.2 


25.4 


353.8 


0.0001572 


0.0010949 


297.5 


0.0000185 





0001980 


HD 190228 b 


1136.1 


0. 


531 


101.2 


40.0 


420.4 


0.0001656 


0.0008165 


353.5 


0.0000274 





0001 936 


CoRoT-10 b 


13.2 


0. 


530 


218.9 


0.0 


1123.2 


0.6402574 


0.9856290 


944.5 


0.3219422 





5607959 


HD 142022 b 


1928.0 


0. 


530 


170.0 


8.3 


273.5 


0.0000047 


0.0000650 


230.0 


0.0000003 


o 


0000(171 


HD 87883 b 


2754.0 


0. 


530 


291.0 


0.0 


193.4 


0.0000000 


0.0000010 


162.6 


0.0000000 





0000000 


HD 108147 b 


10.9 


0. 


530 


308.0 


0.0 


1828.5 


0.7475883 


0.9266767 


1537.6 


0.4740451 





6283575 


HD 168443 b 


58.1 


0. 


529 


172.9 


0.0 


959.0 


0.1119503 


0.1943817 


806.4 


0.0504508 


o 


1 01 fiS^S 

. 1U1UOOO 


HD 81040 b 


1001.7 


0. 


526 


81.3 


20.3 


325.6 


0.0000604 


0.0005207 


273.8 


0.0000059 





.0000813 


HIP 5158 b 


345.7 


0. 


520 


252.0 


49.8 


425.7 


0.0012029 


0.0057508 


358.0 


0.0002036 





.0013882 


HD 148156 b 


1027.0 


0. 


520 


35.0 


24.9 


350.4 


0.0001130 


0.0008256 


294.6 


0.0000131 


0. 


.0001469 


HD 217107 c 


4270.0 


0. 


.517 


198.6 


0.0 


197.0 


0.0000000 


0.0000008 


165.6 


0.0000000 





.0000000 


HAT-P-2 b 


5.6 


0. 


.517 


185.2 


0.0 


2498.6 


1.2828619 


1.4573069 


2101.0 


0.8890871 


1. 


.0564045 


HD 1237 b 


133.7 


0. 


511 


290.7 


42.5 


585.8 


0.0179584 


0.0520765 


492.6 


0.0049349 


0. 


.0184897 


HD 142415 b 


386.3 


0. 


500 


255.0 


72.3 


487.3 


0.0023210 


0.0088971 


409.8 


0.0004917 





.0025685 


HD 215497 c 


567.9 


0. 


490 


45.0 


35.6 


376.9 


0.0002662 


0.0016269 


317.0 


0.0000358 





.0003268 


HD 106252 b 


1531.0 


0. 


482 


292.8 


17.1 


302.4 


0.0000157 


0.0001642 


254.3 


0.0000013 





.0000222 


HD 33636 b 


2127.7 


0. 


481 


339.5 


4.1 


262.2 


0.0000024 


0.0000372 


220.5 


0.0000001 


0. 


.0000037 


HD 181433 d 


2172.0 


0. 


480 


330.0 


0.0 


221.7 


0.0000002 


0.0000057 


186.4 


0.0000000 





.0000004 


HD 196885 b 


1333.0 


0. 


480 


78.0 


35.9 


366.0 


0.0000877 


0.0005804 


307.8 


0.0000111 


0. 


.0001112 


HD 33283 b 


18.2 


0. 


480 


155.8 


0.0 


1503.7 


0.3606475 


0.4817197 


1264.5 


0.2118062 


0. 


.3083335 


HD 210277 b 


442.2 


0. 


476 


119.1 


44.6 


398.9 


0.0005740 


0.0031455 


335.5 


0.0000862 





.0006903 


HD 154857 b 


409.0 


0. 


470 


59.0 


32.1 


575.2 


0.0021955 


0.0065333 


483.7 


0.0005893 


0. 


.0022767 



Note. — The last six columns are for the cases of r\ = and r\ = 1 respectively. 



The top two panels of Figure |3] plot the planetary ef- 
fective temperatures and flux values for all 390 planets 
included in the sample, assuming r\ = 1 and A = 4.5 /zm. 
The planetary temperature plot on the left shows a rea- 
sonably even distribution of temperatures for circular 
orbits, as one may expect. However, notice there is a 
downward trend for eccentricities < 0.5 since these plan- 
ets tend to lie at a large semi-major axis. A second trend 
occurs in the opposite direction (eccentricities > 0.5); the 
temperature increases as the eccentricity increases since 
this decreases the star-planet separation at periastron. 
The temperature plot maps to the flux ratio plot shown 
on the right, where the hot Jupiters can be seen clustered 
in the top-left corner. The upward trend in temperatures 
towards higher eccentricities results in an equivalent up- 



ward trend in flux ratios. There is no decrease in flux 
ratio towards higher eccentricities. Thus it is clear that 
the high eccentricity planets present viable opportuni- 
ties to detect their thermal signature during periastron 
passage. 

The percentage of the total orbital period spent within 
the HZ is calculated by first estimating the boundaries 
of the HZ (see Section [SJ and then determining the star- 
planet separation at equal increments in time during a 
Keplerian orbit. It is not surprising that most of the ec- 
centric planets spend less than half of their time within 
the HZ due to the large range of star-planet separations 
which occur during the orbit. The bottom two panels of 
Figure [3] show how the planets which spend some part of 
their orbit within the HZ arc distributed according to pe- 
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Fig. 3. — Predicted planetary effective temperature at periastron (top-left) and predicted flux ratio at 4.5 fim (top-right) as a function 
of eccentricity. The predicted flux ratios assume a well-mixed (r] = 1) atmospheric model and an albedo of A = 0.0. Also shown are the 
the percentage of the orbit spent within the HZ with respect to orbital period (bottom-left) and the eccentricity as a function of period 
(bottom-right). The last panel only shows planets which enter the HZ where the size of the points linearly increases with the percentage 
of time spent within the HZ. 



riod and eccentricity. The left plot shows the percentage 
of the orbit spent in the HZ, ijjz, as a function of pe- 
riod. Interestingly, the planets which spend more than 
20% of the orbit within the HZ are event distributed 
between 20-80% and orbital periods of 200-2000 days. 
To clarify the relative distributions of the planets which 
spend a non-zero amount of time within the HZ, we show 
the eccentricity of the planets as a function of period in 
the right plot. The relative size of the points increases 
as a function of the percentage time spent within the 
HZ. The planets which spend a portion of their time in 
the HZ are fairly evenly distributed in eccentricity, al- 
though the more circular orbits preferentially spend a 
greater percentage of their time there. This indicates 
eccentricity can sometimes be a useful discriminator in 
selecting targets for potential life-bearing planets, al- 
though there are some relatively large points shown for 
e > 0.5. In both of the bottom plots, the outlier lo- 
cated at 41.4 days is the Saturn-mas s planet orbiting the 
M4 dw arf HIP 57050, discovered by Haghighi pour et all 
(2010). Despite the relatively short period and eccentric- 
ity of 0.31, the smaller luminosity of the host star allows 
this planet to spend most of the orbit within the HZ. 

In Figure S] we show two particularly interesting cases. 
On the left is shown the highly eccentric planet orbiting 



HD 43197. This planet spends ~ 78% of the total orbital 
period residing within the HZ, although even considering 
7/ = 1, the temperature rises to 716 K during periastron 
passage. On the right is shown the rather more benign 
orbit of the planet orbiting HD 156411. This planet is not 
listed in Table [T] (due to its relatively low eccentricity) 
but it spends 96.5% of the orbit within the confines of 
the HZ and reaches a peak temperature of 307 K during 
periastron passage assuming rj = 1. With a flux ratio of 
7.8 x 10 -8 at 4.5 /im, one cannot reasonably expect to 
detect such a planet with current instrumentation, but it 
does present an interesting case for habitability studies 
of eccentric orbits. These two planets perfectly represent 
the two cases described by Figure [5] 

7. IMPACT OF STAR SPOTS ON DETECTION 

Since most of the known exoplanets we are concerned 
with here orbit F, G, or K-type main sequence stars, we 
need to investigate the effect, if any, that star spots will 
have on the detection of thes e planets an d their possible 
phase signatures in the IR. iBerdvuginal ()2005[ ) showed 
that, on average, the difference between the stellar pho- 
tosphere and a stellar spot is larger for hotter stars with 
values near 2000 K for late F and early G stars and drops 
to 200 K for late M stars. IR observations detect stars at 
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Fig. 4. — Habitable zones (shaded region) and orbits of the planets orbiting HD 43197 (left) and HD 156411 (right). 



cooler temperatures than optical observations, and there- 
fore star spots will be observed at a lower contrast than 
in the optical. This means that, even though late-type 
stars typically have a greater number of star spots, these 
should be less of a flux contrast issue for the late-type 
stars that are better suited for IR observations. 

Of concern is the timescale of star spot modulations 
compared with the timescale over which significant flux 
variation is e xpected to occur durin g periastron passage. 
As shown byl.Kane fc Gelinol d20inh . the maximum flux 
amplitude for eccentric planets is a small fraction of the 
total orbit and in many cases allows a 24-48 hour window 
through which the maximum changes can be observed. 
In contrast, the rotation periods of typical exoplanet 
hosting stars is substantially longer. Rotation periods 
h ave been mea s ured f or m any of these, su c h as t he work 
of iHenrv etlll (| 19971 ) and Isimpson et all (|2010t ). These 
rotation periods are mostly in the range of 20-40 days 
which is common of 1-5 Gyr old main sequence stars. 

In the case of the sun, star spots typically cover be- 
tween 10 -3 (during a solar maximum) and 10~ 5 of the 
surface. However, up to a 22% of a hemisphere was seen 
to be covered in a Doppler image of XX Tri, which is a 
K0 giant and there fore not typical of exoplanet host stars 
(|Strassmeierl l2009t ) . Star spots are likely to evolve over 
timescales of a few stellar rotation cycle s or even within 
one cycle. This supports the results of iHussainl (2002) 
who found that spots on single main sequence stars, at 
most, live for weeks. 

The frequency of star spot occurrance increases for low- 
mass stars which, at the current epoch, comprise the 
minority of exoplanet host stars. As described above, 
most stellar rotation periods are much longer than the 
periastron passage timescales. Therefore it is usuaully 
not necessary to worry about star spots inhibiting the 
detection of the phase signatures of the planets. Stel- 
lar spot signatures change over time, so in the very few 
cases where the timescales may be similar, it should be 



possible to disentangle the planet phase signature from 
the star spot signature. C onsistent wit h our fi ndings, 
an IR study conducted by De sert et al.1 (|2011al ) found 
that the variability due to spots is less than the pre- 
dicted t ransit depths and w ith a longer period. Con- 
versely, iDesert et al.1 (|2011b| ) find that they do have to 
take into account spots on HD 189733 and conclude that 
an estimation of the planet-to-star radius ratio should 
be associated with a corresponding stellar activity level. 
Thus, each experiment should address this stellar activ- 
ity level issue on a case- by-case basis. 

8. CASE STUDIES 

In this section we consider several interesting case stud- 
ies from the results in the previous sections. 

8.1. The HD 156846 System 

The planet orbiting HD 156846 has recently been stud- 
ied bv Kane et ail (|2011[ ). providing refined orbital pa- 
rameters and the exclusion of additional giant planets in 
the system. Thus, even though the orbital period is large 
(P = 359 days), observations during periastron passage 
could be timed with great precision in order to detect a 
thermal signature. However, the periastron argument of 
uj = 51° means that the periastron passage occurs close 
to the observer-star line of sight on the near side of the 
star. Thus the night-side of the planet faces towards us 
during periastron passage. This presents an opportunity 
to test the heat redistribution efficiency models for this 
planet since a low efficiency will result in a distinct phase 
function which prevents detection, but a high efficiency 
will still allow for a detectable signature as the heat is 
transferred to the night-side of the planet. 

8.2. The HD 37605 and HD 33283 Systems 

Here we discuss the two interesting cases of HD 37605b 
and HD 33283b. The flux ratio of HD 37605b is second 
only to HD 80606b amongst the eccentric planets shown 



8 Stephen R. Kane & Dawn M. Gelino 




Fig. 5. — Phase curves for HD 37605b (left) and HD 33283b (right) at 4.5 (im, assuming 100% heat redistribution efficiency (r) = 1) of 
the atmosphere. The sub-panels zoom in on the section of the phase curve which would be optimal for monitoring. 



in Table [IJ aided by the short orbital period compared 
with its large eccentricity. This planet was also suggested 
as an excellent c andidate for phase mo nitoring at optical 
wavelengths by iKane fc Gelinol (|2010T) . The periastron 
argument of 10 = 211° means that a significant fraction 
of the day-side will be angled towards us during perias- 
tron passage. HD 33283b has a very similar predicted 
flux ratio, the lower eccentricity being offset by a shorter 
orbital period. The host stars in both cases are bright 
enough such that a detection of phase variations dur- 
ing periastron passage will yield significant information 
regarding the atmospheric properties of the respective 
planets. A lack of detectable phase variation will in turn 
constrain atmospheric models for these planets regarding 
the mechanisms for which heat distribution is occurring. 

Provided the orbital parameters for these stars are suf- 
ficiently refined, the total time needed to achieve ade- 
quate coverage of the phase curves to maximize a detec- 
tion is relatively low. Shown in Figure [5] are the calcu- 
lated flux ratios at 4.5 fim. as a function of orbital phase 
for each of the planets. Zero orbital phase is defined 
to occur at superior conjunction and the location of the 
peak flux ratio depends upon the periastron argument 
of the orbit. This model assumes that the atmospheres 
are very efficient at redistributing the heat r\ ~ 1 which 
lowers the effective flux ratios. Thus, these models rep- 
resent a minimum expected signature from each of the 
planets. Each of the figures contains a sub-panel which 
zooms in on the pericenter passage segment of the orbit. 
These segments can be adequately covered with 25 hours 
of observations, therefore only 50 hours in total would 
be needed to monitor both targets. The planets have 
the same time-coverage requirements because, although 
HD 37605b has a longer period than HD 33283b, it also 
has a larger eccentricity which increases the planetary 
velocity at pericenter. Such targets represent rare oppor- 
tunities to probe the atmospheric properties of planets 
beyond the regime of hot Jupiters and further develop 
theoretical models of these planets. 

9. CONCLUSIONS 

We have presented an analysis of the known exoplan- 
ets to determine predicted temperatures and flux ratios. 
This exercise is clearly not meant to serve as an exhaus- 



tive modeling of the atmospheres fo r each of these planets 
and indeed comb i ning t he models of lFortnev et al.1 ()2010f ) 
and iLewis et al.l ([2010) with IR data will present new 
and exciting opportunities for studying the atmospheres 
of eccentric planets. Rather, this is meant to serve as 
a guide towards which of the particularly eccentric ex- 
oplanets may serve as interesting targets for follow-up 
observations. 

For planets which lie within the HZ of their host 
stars, eccentric planets may present additional opportu- 
nities for studying planetary atmospheres in these zones. 
HD 43197b and HD 156411b are particularly interesting 
cases since they each spend a majority of the total or- 
bital period within their respective habitable zones but 
vary dramatically in the temperature differences experi- 
enced throughout their entire orbits. Thus these cases 
may be used to investigate the effects of eccentric orbits 
upon habitability. Based on the fact that most stellar 
rotation periods are are much longer than the periastron 
passage timescales, and that star spots typically have a 
signature that is much less than that of a habitable zone 
planet at periastron passage, we conclude that spots are 
not a significant issue to this particular study. Due to 
the changing nature of star spot activity on a given star, 
long-term monitoring of the star could disentangle the 
planet phase signature. 

We have presented several case studies from our anal- 
ysis which serve as potentially interesting targets for de- 
tection during periastron passage. HD 156846b is a good 
IR target for studies of the atmospheric properties of 
planets in extreme (highly eccentric) orbits. HD 37605b 
and HD 33283b are also good targets since they have 
bright host stars, relatively high predicted flux ratios, 
and brief periods of large changes in the flux ratio. Warm 
Spitzer is capable of monitoring some of these target dur- 
ing the remaining lifetime of the mission. When this mis- 
sion ceases operations, JWST will have instrumentation 
that is capable of continuing this study to even higher 
precision. 
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